Geophysical and Astrophysical Fluid Dynamics, 1 
Vol. 100, Nos. 4-5, August-October 2006, 341-361 
DOL10. 1 080/03091920600784873 

Parameter dependences of convection-driven 
dynamos in rotating spherical fluid shells 

FRIEDRICH H. BUSSE*f and RADOSTIN D. SIMITEV*| 

* Institute of Physics, University of Bayreuth, D-95440 Bayreuth, Germany 
\ Department of Mathematical Sciences, The University of Liverpool, Liverpool L69 7ZL, UK 

{Received 28 April 2009; in final form ??) 

For the understanding of planetary and stellar dynamos an overview of the major pa- 
rameter dependences of convection driven dynamos in rotating spherical fluid shells 
is desirable. Although the computationally accessible parameter space is limited, ear- 
lier work is extended with emphasis on higher Prandtl numbers and uniform heat 
flux condition at the outer boundary. The transition from dynamos dominated by 
non-axisymmetric components of the magnetic field to those dominated by the ax- 
isymmetric components depends on the magnetic Prandtl number as well as on the 
ordinary Prandtl number for higher values of the rotation parameter r. The depen- 
dence of the transition on the latter parameter is also discussed. A variety of oscillating 
dynamos is presented and interpreted in terms of dynamo waves, standing oscillation 
or modified relaxation oscillations. 

1. Introduction 

Most global magnetic fields of planets and stars are generated by thermal or compo- 
sitional convection in the deep interiors of those bodies. In spite of their common origin 
a large variety of dynamo processes is indicated by the observations of planetary and 
stellar magnetic fields. This is not surprising in view of the widely varying conditions 
under which dynamos operate. It is thus desirable to understand the effects of the most 
influential parameters on convection driven dynamos. The purpose of a study of the 
parameter dependence of dynamos would be twofold: On the one hand properties of ob- 
served magnetic fields could be explained in terms of parameter values of the respective 
system. On the other hand, unknown conditions in the interior of the celestial bodies 
may be inferred from the spatio-temporal structures of their magnetic fields. While it 
is not yet possible to create convincingly detailed models of planetary dynamos and of 
the solar cycle, some major variations of dynamos as function of their parameters can be 
explored through numerical simulations. 

The increasing availability in recent years of computer capacity has facilitated large 
scale numerical simulations of the generation of magnetic fields by convection in rotating 
spherical fluid shells. Because of limited numerical resolution, molecular values of mate- 
rial properties are usually not attainable in computer simulations and eddy diffusivities 
representing the effects of the unresolved scales of the turbulent velocity field must there- 
for be invoked for comparisons with observations. It is often assumed for this reason that 
the eddy diffusivities for velocities, temperature and magnetic fields are identical. The 
effects of turbulence on the diffusion of vector and scalar quantities differ, however, and 
eddy diffusivity ratios such as the effective Prandtl number and the effective magnetic 
Prandtl number thus do not equal unity in general. Besides the Prandtl numbers, the 
boundary conditions exert a strong influence on convection and its dynamo action. Both 
effects will be considered in this paper. 
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In some respects this paper represents an extension of an earlier paper (Simitev and 
Busse, 2005, to which we shall refer to by SB05) which has focused on the dependence of 
average properties of convection driven dynamos on the Prandtl number. New results will 
be reported in the following and the emphasis will be placed on time dependent properties 
and on dynamo oscillations in particular. In addition the effect of boundaries of low 
thermal conductivity will be studied which arc often more realistic than the commonly 
assumed boundaries with fixed temperature. Especially in the case of the Earth the low 
conductivity of the mantle will lead to a uniform heat flux from the core unless the effect 
of mantle convection is taken into account. Since the inhomogeneity introduced by the 
latter is not well known it will not be considered in the present analysis. For simulations 
with various inhomogeneous thermal boundary conditions see the papers by Glatzmaicr 
et al. (1999) and by Olson and Christensen (2002). If compositionally driven convection 
in the Earth's core is emphasized the choice of uniform flux at the inner boundary and 
of fixed composition at the outer boundary is appropriate as has been assumed in the 
work of Glatzmaicr and Roberts (1995). 

After a brief introduction of the basic equations and the method of their numerical 
solution in section 2, some properties of convection without magnetic field will be consid- 
ered in section 3. In particular the influence of fixed heat flux boundary conditions will 
be explored. In section 4 the onset of convection driven dynamos in fluids with different 
Prandtl numbers is described and in section 5 the oscillatory dynamos are interpreted in 
terms of the theory of dynamo waves. The influences of various boundary conditions are 
studied in section 6 and a concluding discussion is given in the final section 7. 

2. Mathematical formulation of the problem and methods of solution 

We consider a rotating spherical fluid shell of thickness d and assume that a static 
state exists with the temperature distribution T$ = To — f3d 2 r 2 /2. Here (3 = q/(3Kc p ) 
and T = T\ — (T 2 — Ti)/(1 — rf), where T\ and T 2 are the constant temperatures at 
the inner and outer spherical boundaries, rj = r,/r is the radius ratio of the inner r, 
to the outer r a radius, q is the uniform heat source density, k its thermal diffusivity, c p 
is its specific heat at constant pressure and rd is the length of the position vector with 
respect to the center of the sphere. The gravity field is given by g — —djr. In addition 
to d, the time d 2 /v, the temperature v 2 /jad 4 and the magnetic flux density u(fig) 1 ^ 2 jd 
are used as scales for the dimcnsionlcss description of the problem where v denotes the 
kinematic viscosity of the fluid, g its density and fi is its magnetic permeability. Since we 
shall assume the Boussinesq approximation material properties are regarded as constants 
except for the temperature dependence of the density described by a = — (dg/ dT)/g 
which is taken in the gravity term. Both, the velocity field u and the magnetic flux 
density B, arc solenoidal vector fields for which the general representation 



can be employed. By multiplying the (curl) 2 and the curl of the Navier-Stokes equations 
of motion by r wc obtain two equations for v and w, 



[(V 2 - d t )C 2 + rd v }X7 2 v + tQw - £ 2 6 = -r ■ V x [V x (u ■ Vu - B ■ VB)], (2.2a) 



where dt denotes the partial derivative with respect to time t and where d v is the 



u = V x (Vv x r) + Vro x r, 
B = Vx (Vft x r) + Vg x r, 



(2.1a) 
(2.16) 



[(V 2 - d t )C 2 + rd v ]w - tQv = r ■ V x (it ■ Vu - B ■ VB) 



(2.26) 
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Figure 1. Critical values of the Rayleigh number R c , wavenumber m c and frequency ui c as 
functions of r in the case n = 0.4 for values of the Prandtl number P = 0.2 (red), P — 2 (blue), 
P = 20 (green) and UFBC (thick lines, filled symbols) and FTBC (thin lines, empty symbols). 



partial derivative with respect to the angle tp of a spherical system of coordinates r, 9, ip. 
For further details we refer to SB05. The operators C 2 and Q are defined by 

£ 2 = -r 2 V 2 + d r (r 2 d r ), 

Q = r cos(9V 2 - (£ 2 + rd r )(cos6d r - r' 1 sm0dg). 

The heat equation for the dimensionlcss deviation from the static temperature distri- 
bution can be written in the form 

V 2 6 + RC 2 v = P{d t + u ■ V)6, (2.3) 

and the equations for h and g are obtained through the multiplication of equation of 
induction and of its curl by r 

V 2 C 2 h = P m [d t C 2 h -fVx(uxB)], (2.4a) 
V 2 £ 2 <? = P m [d t C 2 g - r • V x (V x (it x B))}. (2.46) 

The Rayleigh number R, the Coriolis number r, the Prandtl number P and the magnetic 
Prandtl number P m are defined by 

^ K A 

where A is the magnetic diffusivity. For the static temperature distribution we have 
chosen the case of a homogeneously heated sphere. This state is traditionally used for 
the analysis of convection in self-gravitating spheres and offers the numerical advantage 
that for Rayleigh numbers close to the critical value R c the strength of convection does 
not differ much near the inner and outer boundaries. As the Rayleigh number increases 
beyond R c additional heat enters at the inner boundary and is delivered by convection to 
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Figure 2. Time series of kinetic E energy densities (top) and Nusselt numbers Nu (bottom) 
for the convection case P = 0.5, t = 1.5 x 10 4 , R = 1.5 x 10 6 and UFBC. The components E± 
and Nu are represented by solid, Et and Nut by dotted lines and E v by a dashed line. 




Figure 3. Isolines of azimuthally-averaged differential rotation u v (left half) and temperature 
perturbation 9 (right half) in the case P = 20, r = 10 4 , R = 2 x 10 6 and FTBC (to the left) or 
UFBC (to the right). The fields represent averages in time. 



the outer boundary. When R reaches a high multiple of R c the heat generated internally 
in the fluid becomes negligible in comparison to the heat transported by convection 
through the spherical shell. 

For the velocity field either stress-free boundaries, 

v = d rr v = d r (w/r) = at r = j-j = 77/(1 — 77) and r = r a = 1/(1 — 77), (2.6) 

or no-slip boundaries, 

v = d r v = w = at r = ^ = 77/(1 — 77) and r = r = 1/(1 — 77), (2-7) 

will be assumed. Unless indicated otherwise 77 = 0.4 and stress-free boundary conditions 
will be used in the following. The value 77 = 0.4 is not only employed because it turns 
out to be close to the radius ratio of the liquid outer core of the Earth, but because it 
also provides a good balance between the regions inside and outside the tangent cylinder 
touching the inner boundary at its equator. These two regions are known, of course, for 
their rather different dynamical properties. 

For the thermal boundary conditions we shall assume that either the temperature is 
fixed (FTBC, i.e. Fixed Temperature Boundary Conditions), 

9 = (2.8) 



or its normal derivative is fixed (UFBC, i.e. Uniform Flux Boundary Conditions), 
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Figure 4. Time- and azimuthally-averaged temperature 0(r ,9) in the case UFBC (left plot 
(a)) and Nusselt number Nu (9) in the case FTBC (right plot (b)) at the outer boundary as a 
function of the colatitude 6 for P = 0.5 with r = 1.5 x 10 4 and R = (10 + 5ra) x 10 5 ,n = 0...4 from 
bottom to top at 0° (solid lines) and for P = 20 with r = 10 4 and R = (10 + 5n) x 10 5 , n = 0...6 
from bottom to top at 90° (dashed lines). The temperature values Q(r o ,0) for P = 0.5 have 
been divided by 20. 



at the boundary where the double overbar indicates the average over the spherical bound- 
ary. This latter definition has been used in order that the Rayleigh number continues to 
be based on the difference between the average temperatures between the boundaries. 

For the magnetic field electrically insulating boundaries are assumed such that the 
poloidal function h must be matched to the function which describes the potential 
fields outside the fluid shell 

g = h-h {e) = d r (h-h {e) ) = Q at r = 77/(1 - rj) and r = 1/(1 -77). (2.10) 

But computations for the case of an inner boundary with no-slip conditions and an 
electrical conductivity equal to that of the fluid have also been done. The numerical 
integration of equations ()2.2j) . (|2.3j) and (|2.4p together with boundary conditions (|2.6|) or 
(|2.7|) and (|2.8|) or (|2.9p and (|2.10p proceeds with the pseudo-spectral method as described 
by Tilgner and Busse (1997) and Tilgner (1999) which is based on an expansion of all 
dependent variables in spherical harmonics for the 9, (^-dependences, i.e. 

v = Vi" 1 {r,t)P™ {cos 8) exp{imp} (2.11) 

Lin 

and analogous expressions for the other variables, w, 0, h and g. P™ denotes the asso- 
ciated Lcgcndrc functions. For the r-dcpcndcncc expansions in Chcbychcv polynomials 
are used. For further details see also Busse et al. (1998) or Grote et al. (2000). For the 
computations to be reported in the following a minimum of 33 collocation points in the 
radial direction and spherical harmonics up to the order 64 have been used. But in many 
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Figure 5. Kinetic E energy densities (left ordinate) and Nu (right ordinate) as functions of 
R for non-magnetic convection for (a) P = 0.5, r = 1.5 x 10 4 and (b) P = 20, r = 10 4 . The 
components Et, E p , Et are represented by squares, plus-signs and crosses, respectively and Nu 
and Nm by diamonds and triangles. Cases with FTBC are shown with light symbols, and these 
with UFBC are shown with heavy symbols. 



cases the resolution has been increased to 49 collocation points and spherical harmonics 
up to the order 96 or 128. 



3. Convection in rotating spherical shells 

For an introduction to the problem of convection in spherical shells we refer to the 
recent review of Busse (2002a). Additional information can be found in the papers by 
Grote and Busse (2001), Christcnscn (2002) and Simitev and Busse (2003, 2005). Here we 
shall focus the attention on the case of uniform heat flux boundary conditions. In figure 
[1] the critical value R c of the Rayleigh number and the associated critical azimuthal 
wavenumbcr m c have been plotted as a function of r for several values of the Prandtl 
number P. The results for the case of fixed temperatures are shown for comparison. In 
the case of the steady onset of convection in a planar layer heated from below the critical 
value of the Rayleigh number and the associated wavenumber decrease monotonously 
with a decreasing thermal conductivity of a boundary (Sparrow et ai, 1964). There the 
onset of convection is described by a one-dimensional eigenvalue problem, while in the 
present case a two-dimensional problem in the complex domain must be solved because of 
the time dependent onset. Hence exceptions from the monotonous decrease of the critical 
values R c , m c in going from the fixed temperature to fixed flux boundary condition must 
be expected as are indeed found in the intermediate range of r, 10 3 < r < 3 x 10 4 . This 
phenomenon can also be recognized in the analytical results of Busse and Simitev (2004) 
obtained in the low Prandtl number limit. As the critical wavenumber m c increases the 
value R c becomes independent of the boundary conditions as must be anticipated since 
the r-dependence becomes negligible in comparison with the azimuthal dependence of 
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Figure 6. Kinetic E x and magnetic M x energy densities (left) and viscous V x and Ohmic O x 
dissipations (right) as functions of P for convection driven dynamos for r = 10 4 , R = 2 x 10 6 , 
FTBC and magnetic Prandtl number as indicated in the boxes. The components X p , X t , X p , Xt 
(where X = E, M, V , O) are represented by circles, squares, plus-signs and crosses, respectively. 
Kinetic energy densities and viscous dissipations are shown with light symbols, magnetic energy 
densities and Ohmic dissipations are shown with heavy symbols. E p is multiplied by a factor 
20. 



convection. The same independence is, of course, also evident in the asymptotic theory 
of Bussc (1970). 

At finite amplitudes convection in the presence of uniform heat flux boundaries evolves 
in a similar way as convection with fixed temperatures at the boundaries as long as the 
wavenumbers m c of the two cases do not differ significantly. It should be mentioned here 
that in contrast to the linear results displayed in figure Q] only an outer uniform-flux 
boundary has been assumed for the computations of finite amplitude convection and its 
dynamos, while the fixed temperature condition has been kept at the inner boundary. This 
combination seems to be most appropriate for some planetary applications as mentioned 
in the Introduction. 

The initial regime of drifting thermal Rossby waves is followed by the subsequent 
stages of vacillating convection, localized convection and relaxation oscillation that are 
observed with increasing Rayleigh number for Prandtl numbers of the order unity or less 
(but above values of P where convection in the form of inertial waves dominates, i.e. 
P > 10/y/r according to Ardes et al. (1996)). No significant influence of the thermal 
boundary conditions can be noticed. The relaxation oscillations shown in figure [2] seem 
to be a little more pronounced than in the case of a fixed temperature outer boundary, 
but are otherwise identical. At higher Prandtl numbers a stronger tendency towards the 
thermal wind type of differential rotation can be noticed in the presence of a uniform-flux 
thermal boundary as seen in figure [3l This is a consequence of the stronger latitudinal 
temperature gradient admitted by the insulating boundary condition. 

At finite amplitudes of convection the heat transport is of special interest the efficiency 
of which is measured by the Nusselt number. This is defined as the heat transport in 
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Figure 7. Convection-driven dynamos as function of R, P and P m for r = 5 x 10 3 . The symbols 
indicate dynamos for which M p < M p (V), those for which opposite is true (A) and decaying 
dynamos (0). 



the presence of convection divided by the heat transport in the absence of motion. In 
the case of the spherical fluid shell two Nusselt numbers can be defined measuring the 
efficiency of convection at the inner and the outer boundary, 



1 



p de 

Ti dr 



Nu n 



1 



P dB 

r a dr 



(3.1) 



where the double bar indicates the average over the spherical surface. In addition the 
local Nusselt numbers as function of latitude 



Nui{0) = 1 - 



p de 

Tj dr 



Nu (9) = 1 - 



P dB 

r„ dr 



(3.2) 



are of interest where only the azimuthal average is applied, as indicated by the single 
bar. The latitudinal dependence of the local Nusselt number and of the azimuthally 
averaged temperature at the outer boundary are shown in figure 0] for the cases of fixed 
temperature conditions (right) and uniform flux conditions (left) in the cases P = 0.5 and 
P = 20. As must be expected latitudes of high temperatures at the boundary correspond 
to those with high values of Nu o (0). This figure demonstrates that at low supercritical 
Rayleigh numbers the heat transport occurs primarily across the equatorial region, but as 
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Figure 8. Convection-driven dynamos as function of P and P m with r = 10 4 and R — 2 x 10 6 
(light symbols) and with t — 3 x 10 4 and R = 3.5 x 10 6 (heavy symbols). The dynamos for 
which M p > M p are indicated by squares in the case of r = 3 x 10 4 and by plus-signs in the 
case of r = 10 4 . The dynamos for which M p < M„ are indicated by diamonds in the case of 
r = 3 x 10 4 and by crosses in the case of r = 10 . The broken and the solid lines represent 
the approximate location of the parameter values where the transition from M p -dominated to 
Mp-dominated dynamos occurs for t — 3 x 10 4 and r = 10 4 , respectively. 



R increases the heat transport in the polar regions takes off, especially in the case of the 
lower value of P, and soon exceeds that at low latitudes. In the polar regions convection 
is better adjusted for carrying heat from the lower boundary to the upper one, and it is 
known from computations of the convective heat transport in horizontal layers rotating 
about a vertical axis that the value of Nu may exceed the value in a non-rotating layer 
at a given value of R in spite of the higher critical Rayleigh number in the former case. 

At higher values of the Prandtl number the takeover of the heat transport by convection 
in the polar regions is less pronounced, but the contrast between the equatorial region and 
the mid-latitudes is even stronger than for lower values of P. Again, a dramatic difference 
between fixed flux and fixed temperature boundary conditions can not be discerned. 

It is also useful to compare the averages over space and time of the kinetic energy 
densities of the various components of the convection flow. They are defined by 

E P = i(| V x (V« x r) | 2 ), E t = ^(1 Vwxr | 2 ), (3.3a) 

4 = J(|Vx(Vi;xr) | 2 >, E t = i(| W; x r | 2 ), (3.36) 

where the angular brackets indicate the average over the fluid shell and over time and v 
refers to the azimuthally averaged component of v, while v is defined by v = v — v. In 
figure [5] the growth of these energies with the the Rayleigh number is shown for both, 
the case with an insulating and the case with a fixed temperature at the outer boundary. 
In the low Prandtl number case of the left diagram the rapid growth with R of E t 
corresponding to the energy of differential rotation is remarkable. Only for higher values 
of P does Et never exceed the energies of the fluctuating components of motion. 
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Figure 9. Time-averaged kinetic energy densities E p , (solid line and circles), Et (dashed line 
and squares), E p (dashed-dotted line and plus-signs), Et (dotted line and crosses) and Nusselt 
number Nui (thin line and solid diamonds) of non-magnetic convection (lines) and dynamo 
solutions (symbols) for P — 20, r = 5 x 10 3 and R as indicated at the x-axis. The symbols 
connected by the dotted lines correspond to the values of P m given at the top of the figure. 

4. Convection-driven dynamos 

The onset convection-driven dynamos as function of parameters equivalent to those 
given by expressions (2.5) has been investigated in numerous papers, see, for example, 
Busse et al. (1998), Christensen et al. (1999), Grote et al. (2000, 2001), Grote and Busse 
(2001), Kutzner and Christensen (2000, 2002), SB05. For displays of dynamo onsets 
as functions of P,R and r see SB05 and Busse (2002). Typically the Raylcigh number 
required for dynamos increases with decreasing P m such that the minimum magnetic 
Reynolds number Rm defined by Rm = (2E) 1 / 2 P m stays constant at a value of the 
order 10 2 . As has been found before (Christensen et al. 1999; SB05) a too high value 
of Rm can be detrimental for dynamos because of flux expulsion from the convection 
eddies. Especially at Prandtl numbers of the order unity or somewhat lower a finite 
amplitude of the magnetic field is required for dynamos close to their minimum value of 
R since the velocity field modified by the action of the Lorentz force is more suitable for 
dynamo action, mainly through the suppression of the differential rotation (Grote and 
Busse, 2001), than the non-magnetic velocity field. In a statistical sense convection-driven 
dynamos thus usually correspond to saddle node bifurcations. 

In SB05 it has been noted that a transition from dynamos dominated by non-axisymmetric 
components of the magnetic field to those dominated by the axisymmetric poloidal field 
occurs with increasing Prandtl number when a value P = Pt ~ 8 is reached in the case 
t = 5 x 10 3 . In figure [5] results are displayed for the case r = 10 4 which also exhibit this 
transition. The magnetic energies M p , M t , M p , Mt are defined in analogy to the defini- 
tions (3) for the kinetic energies with h and g replacing v and w. For P < P t , M p < M p 
holds and for P > P t the opposite is true. The results of figure [6] demonstrate that the 
transition Prandtl number Pt depends strongly on P m with higher values of P m favoring 
higher values of Pt- For this reason a more complete analysis of the transition in the 
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Figure 10. A period of dipolar oscillations for P = 5, r = 5 x 10 3 , R = 6 x 10 5 , Pm = 10 with 
A = 0.1429. The plots start at the upper left corner and follow clockwise with At = 0.035 and 
show meridional lines of constant B v in the left half and of rsmOdgh =const. in the right half 
of each circle. 

case t = 5x 10 3 has been made with the inclusion of the Rayleigh number. The results 
shown in figure [7] indicate that the parameter P m appears to have a lesser influence on 
Pt at r = 5 x 10 3 than at r = 10 4 . The trend that with increasing r the value of Pt 
decreases is confirmed in figure [8] where results obtained for r = 3 x 10 4 and for r = 10 4 
are compared. In each of the two cases the Rayleigh number R has been kept constant. 
While the energies of the axisymmctric components of the magnetic field increase relative 
to those of the non-axisymmetric components with increasing P, only a weak tendency 
in this direction can be noticed in the corresponding Ohmic dissipations. 

In SB05 the increase in the efficiency of convection in the presence of a dynamo gen- 
erated magnetic field has been demonstrated for Prandtl numbers of the order unity or 
less. At higher Prandtl numbers the magnetic field is no longer needed to oppose the 
destructive (for the convective heat transport) effects of the differential rotation which 
manifest themselves, for instance, in the relaxation oscillations shown in figure [5J When 
Nusselt numbers and kinetic energies are thus compared for P = 20 in figurc[5]for convec- 
tion with and without a dynamo, little change is seen. The energy of differential rotation 
even appears to be slightly larger in the magnetic case. 

5. Oscillatory dynamos 

It has long been known that mean-field dynamos in the presence of a strong differential 
rotation typically occur in the form of oscillatory dynamos, especially if any meridional 
circulation is relatively weak (Roberts, 1972). In the present case of dynamos driven 
by turbulent convection nearly time periodic oscillations are a persistent feature in the 
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Figure 11. Period of dynamo oscillations T obtained from the analytical expression (5.8) 
(empty symbols) and from simulations (solid symbols) for (a) P = 0.1, r = 10 s , (b) P = 1, 
t = 10 4 , (c) P = 1, r = 3 x 10 4 and (d) P = 1, (e) P = 3, (f) P = 5 and r = 5 x 10 3 . The values 
of R are given at the x-axis and the values of P m are given at the top of the panels. The symbols 
indicate quadrupolar (circles), hemispherical (squares) and dipolar (triangles) dynamos. 



case of hemispherical and quadrupolar dynamos. Typical examples have been presented 
by Grote et al. (1999) and Grote and Busse (2000). Oscillatory dynamos with dipolar 
symmetry are less frequently encountered. An example is shown in figure 1101 Since the 
prograde differential rotation increases with radius r the oscillations assume the form of 
two symmetric dynamo waves propagating from the equator to high latitudes where they 
decay similarly as found in mean-field dynamos. 

For the theoretical description of these waves we follow Parker's (1955) plane layer 
model with the (x, y, z)-coordinates corresponding to the (tp, — 9, redirections in the 
spherical shell. With the ansatz 

B = B p + iB, B p = VxiA, v = iU + v, (5.1) 

where i is the unit vector in the x-direction and where the fields A, B are rr-indepcndcnt, 
U depends only on z and the x, y-average of v vanishes. Using Parker's equations in 
dimcnsionless form, 

-A = aB + V 2 A/P m , -B = B p -VU + V 2 B/P m , (5.2) 
we may obtain solutions of the form 



(A, B) = (A, B) exp[iq ■ x + at}. 



(5.3) 
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The equations are satisfied with 

pA = aB, pB = -i(q x VU) X A with p = a + \q\ 2 /P m (5.4) 
and with the resulting dispersion relation 

p 2 = 2iT = 2i(-a(q x VU) x /2). (5.5) 
The latter yields the growth rates 

f -\ q \ 2 /p rn ± Vr± iVT for r>o, 

a = < , , (5.6) 

\ -| 9 | 2 /p m ±^TlTi^r| for r<o. 

Growing waves are found for =paq y dU/dz > 2\q\ 4 /P% l with the wave propagating in the 
=py-dircction provided that q y and dU/dz are positive quantities. Here the upper (lower) 
sign refers to the case T > (r < 0). 

The coefficient a is related to the helicity of the convection flow. 

(5.7) 

where F{q,u>) is the helicity spectrum function (Moffatt. 1978) which is assumed to have 
a peak at the typical wavenumber q with q 2 /P m \lu\. q will typically exceed the value 
\q\ used in expression (5.3). But if we use an average value q y ps q ps 2n/d, we find for 
the dimcnsionlcss frequency uj, 

1 / 7T 



P m -(v-V xv)\/2E t ) . (5. 
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Figure 13. Same as figurefTDIbut for P = 0.05, r = 10 5 , R = 4x 10 6 , P m = 0.3 with A = 1.1364. 
with At = 0.0035. The time series of plots starts with the upper left plot and continues in the 
clockwise sense. 

Periods corresponding to this value of u> are presented for several dipolar, hemispherical 
and quadrupolar dynamos in figure [IT] where they can be compared with periods obtained 
from the dynamo simulations. For this purpose the helicity (v ■ V x v) was obtained as an 
average over time and over either the northern or the southern hemisphere. Even though 
the expression (5.8) represents only a crude estimate it seems to fit the various oscillatory 
dynamos quite well. The agreement between mean-field theory and numerical simulation 
is due to the dominant action of the differential rotation in creating the axisymmctric 
toroidal magnetic flux. As discussed in more detail in SB05 the a-effect of mean-field 
theory is usually not clearly noticeable in numerical simulations in contrast to the uj- 
effect (see also Schrinner et at. 2005) 

Because the Prandtl number is moderately large in the case of figure [TO] the typical 
polar flux tubes have formed which do not participate in the oscillation. The axisymmetric 
poloidal field, however, participates fully in the oscillation and exhibits about equal 
amplitudes for both polarities. This property is often missing as shown in figure [121 where 
the red field lines dominate and a dipole field with the opposite polarity does not become 
established throughout the cycle. Finally, an extreme version of the latter situation is 
found in the case of the so-called "invisibly" oscillating dynamos (Busse and Simitev, 
2005). A typical example is found when the magnetic Prandtl number is somewhat lower 
than in the case of figure [TO] such that the polar flux tubes are much stronger and the 
oscillation of the axisymmctric poloidal field is suppressed. As a result external observers 
will not be able to recognize the oscillation of the axisymmetric toroidal field outside the 
tangent cylinder. 

Besides the oscillations related to traveling dynamo waves other types of oscillations 
may be noticed. Typical for Prandtl numbers of the order unity or less are global os- 
cillations in which a nearly hemispherical field switches periodically the hemispheres as 
shown in figure 1131 This behavior reminds one of similar, but much weaker oscillations 
with a period of about 3.6 years which have been reported for the solar magnetic field 
(Knaack et al. 2004). According to the observations the flux emerging from the solar 
southern hemisphere periodically exceeds that from the northern hemisphere and vice 
versa. 

While the dynamo oscillations can be well understood in terms of the linear Parker 
model, they actually occur at finite amplitudes of the magnetic field. As a measure of 
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Figure 14. A period of hemispherical oscillations for P = 0.1, r = 10 5 , R = 6 x 10 6 , P m = 0.11 
with A = 0.1424. The plots follow clockwise with At = 0.005 and show isolines similar to those 
in figure ITuI 




0.05 0.07 0.09 0.11 0.13 



Figure 15. Time series of energy densities of a hemispherical dynamo in the case P = 0.1, 
r = 10 5 , R = 6 x 10 6 , P m = 0.11. The upper and middle panels show energy densities of dipolar 
and quadrupolar components of the magnetic field, while the lower panel shows energy densities 
of the velocity field. The mean toroidal components are represented by solid lines, the fluctuating 
toroidal - by dotted lines, the mean poloidal - by dot-dashed lines and the fluctuating poloidal 
by dashed lines. 



the latter we have denoted the Elsasser number A in the captions of figures QUI 021 H31 
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Figure 16. The dynamo cases (a) to (h) presented in table [T] (from left to right, first row then 
second row). In each plot the upper left quarter shows the differential rotation u v , the lower left 
quarter shows the temperature perturbation O and the right half shows the meridional field lines 
rsmOdgh =const. The plot of case (f) shows the meridional streamlines rsinOdgv =const. to 
the right. 



and [14] The Elsasser number is defined by 

2MP m , s 

A = (5.9) 

T 

where M denotes the magnetic energy density averaged over the fluid shell and in time, 
i.e. M = M p + Mt + M p + M t . The interaction between magnetic and velocity field has 
been taken into account in expression (5.8), of course, in that the actual values for E t 
and the helicity have been used in its evaluation. An even stronger interaction with the 
velocity field is found in the case of figure [14] where a coupling with relaxation oscillations 
of convection similar to those shown in figure [2] occurs. The period, however, is still 
comparable to that of the dynamo waves discussed above. These oscillations are typical 
for hemispherical dynamos of the kind shown in figure [T4l as is evident from figure [TlT a). 
Hardly a full wavelength of the standing dynamo wave participates in the oscillation in 
this case and thus the amplitude of the mean poloidal field varies throughout the cycle. 
Accordingly the Lorentz force restraining the differential rotation varies as well. The 
differential rotation thus participates in the oscillation and the other components of the 
velocity field affected by its shearing action do so as well as can clearly be seen in the 
time record of the kinetic and magnetic energies displayed in figure [15] for the same case 
as shown in figure [T4l 



6. Dynamos for different boundary conditions 

All dynamos discussed so far in the preceding sections have been computed for stress- 
free boundaries with fixed temperatures. In this section some convection-driven dynamos 
obtained for other boundary conditions are considered. They have been listed in table 
[T] As must be expected on the basis of the results shown in figure [5] the strength of 
convection with fixed heat flux outer boundary is lower than in the case of the fixed 
temperature condition. This property is reflected in the corresponding dynamos where 
the kinetic and magnetic energies in cases B arc usually less than those in cases A except 
for the differential rotation and the associated mean toroidal magnetic field. 

Of particular interest are the comparisons of cases with and without magnetic fields. 
Values for the latter cases have been added in table 1 in brackets wherever they were 
available. In the cases (a) and (b) the numerical resolution turned out to be insufficient in 
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Table 1 . Time-averaged global properties of dynamos with various velocity and thermal bound- 
ary conditions as follows. A: stress-free and fixed-temperature, B: stress-free and fixed-flux, C: 
no-slip and fixed-temperature, D: no-slip and fixed- flux. The predominant symmetry type is 
indicated with "D" if dipolar, "Q" if quadrupolar and "— " if the dynamo is decaying. Values for 
non-magnetic convection are given in brackets where available. 



the absence of the magnetic field such that the solution diverges. This result emphasizes 
the smoothing effect of the magnetic field in cases with P m < 1. At high Prandtl numbers 
the kinetic energies in dynamo cases are always a bit lower than in the corresponding 
non-magnetic cases. Only the mean poloidal energy densities measuring the strength 
of axisymmctric meridional circulations and in some instances the mean toroidal energy 
density appear to be enhanced by the dynamo action. This result seems to be independent 
of the thermal boundary conditions and agrees generally with the results displayed in 
figure [9l This same situation appears to hold in the case of no-slip conditions at the 
boundaries even for Prandtl numbers of the order unity or less as is evident from case 
(c) of table 1. This property is in stark contrast to the case of stress-free boundaries 
where the energy densities of the non-axisymmetric components are usually strongly 
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amplified by the dynamo as is evident from figure 16 and others of SB05. Only the 
energy of differential rotation is reduced by the Lorentz force in both, the no-slip and 
the stress-free cases. The same observations apply in the comparison of results obtained 
by Christensen et al. (1999) for the two types of boundary conditions. 

The analogy between the results for no-slip boundaries and for stress-free boundaries 
at high Prandtl numbers is not surprising if it is remembered that the presence of Ek- 
man layers at the rigid boundaries increases strongly the viscous dissipation and thus 
the effective Prandtl number relative to the stress-free case. The analogy holds also for 
the energy of the axisymmetric dipolar field which exceeds the energies of all other com- 
ponents of the magnetic field in the no-slip cases (c) and (d) as well as in the stress-free 
cases (g) and (f) at P = 20 in table 1. 

Characteristic structures of velocity and magnetic fields of the examples of table 1 
are shown in figure 1161 Only a quarter of the respective circles is used to visualize the 
differential rotation and the axisymmetric component of O since these fields are nearly 
symmetric relative to the equatorial plane. 

7. Concluding remarks 

It is popular in dynamo theory to neglect the momentum advection terms in the equa- 
tion of motion (Glatzmaier and Roberts, 1995a; Jones, 2003) since they are undoubtedly 
small in comparison with the Coriolis force. In view of the results of section 4 and the 
analysis of SB05 this neglect can be well justified in the limit of high Prandtl numbers and 
to some extent in the presence of strong influences of Ekman layers. It is not surprising 
that dynamos without momentum advection terms exhibit the typical features of high 
Prandtl number dynamos such as a dominant axial dipolar magnetic field with strong 
toroidal flux tubes in the polar regions. But for values of P of the order unity or lower 
and stress-free boundary conditions the neglect of the momentum advection terms is less 
well satisfied since quite different types of dynamos with dominant non-axisymmetric 
components of the field are found as seen in figure [6] and even more dramatically in cor- 
responding figures of SB05. The transition to dynamos with dominating higher multipoles 
found by Kutzner and Christensen (2002) at higher Rayleigh numbers in the presence of 
no-slip boundaries must probably be attributed also to the momentum advection terms. 
In spite of their smallncss they can drive the geostrophic differential rotation which is the 
most easily exited component of fluid motion. Neither Coriolis force, nor pressure gra- 
dient, nor buoyancy force can drive a geostrophic differential rotation. Only the Lorentz 
force enters as a possible competitor. The latter force usually tends to damp the differ- 
ential rotation according to Lcnz' rule since the mean toroidal field is generated by it. It 
seems that the separation of dynamos into these two classes will persist when r increases. 
But further investigations of this point are desirable. 

The analysis of this paper has focused on cases of moderate Rayleigh numbers where 
the effects of the particular choice of the basic temperature profile are still noticeable. 
Other authors have employed temperature profiles without internal heating and Kutzner 
and Christensen (2000) have analyzed ways in which different basic profiles affect prop- 
erties of convection driven dynamos. Similarly a fixed radius ratio has been used in the 
present paper in order to avoid the discussion of even more parameters. While convection 
driven dynamos in rotating spherical shells with a radius ratio of 0.4 or lower typically 
exhibit dipole and quadrupolc components of the magnetic field that arc axially oriented, 
this situation tends to change as higher values of the radius ratio are used. Depending on 
initial conditions Aubcrt and Wicht (2004) find in the case of a radius ratio 0.5 dynamos 
with an equatorial dipole or with an axial dipole for identical values of the parameters 
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of the problem. Clearly, in applications to cases of particular planetary or stellar mag- 
netic fields parameters that have been fixed in the present analysis are likely to play an 
important role. 
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